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Abstract 

Stochastic quantization is used to derive exact equations, connecting multilocal 
field correlators in the ip 3 theory and gluodynanrics. Perturbative expansion of the 
obtained equations in the lowest orders is presented. 


1 Introduction 

Exact determination of nonperturbative effects is one of the main problems in field theory 
nowadays. We suggest to look for them in exact equations (eqs.) for vacuum correlators, 
which we derive, using stochastic quantization method [1]. 

In QCD the nonperturbative correlators are an essential part of operator product 
expansion (OPE), QCD sum rules [2] and, more recently, of the method of vacuum corre¬ 
lators (MVC) [3]. In the latter correlators are used as a dynamical input, and one needs 
to obtain eqs. for them, starting from the lagrangian. 

There exist in literature two types of exact eqs. First is Dyson-Schwinger eqs. [4], 
which, however, are not applicable to QCD in the confining phase, because of the lack of 
gauge invariance. The other one is Makeenko-Migdal eqs.[5], which are closer to our eqs., 
but written not for correlators, but for Wilson loops in the large N c limit. We hope, that 
the eqs., found in this letter for correlators, can be used more directly for the purposes, 
stated above. 

Sections 2 and 3 are devoted to the ip 3 theory, and in Section 4 an alternative approach 
is applied to gluodynamics. 

2 (£> 3 theory 

From the Euclidean action 

S = J dx 
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in the stochastic quantization method one obtains the Langevin equation 

(p(x, t) + ( m 2 - d 2 - g(p(x, t ) = rj(x, t) , (1) 

where < g{x, t)rj(x', t') > = 2 5(x — x')5(t — t'). In the absence of constant classical 

solutions the retarded Green function of (1) is 

G(x,y,t) = 9{t) < x\e- (m2 - a2 ~ 9(p)t \y > = 


where 



( Dz) xy — 


N 


lim 

N— KX) 


n 

n =1 


d A z{n) 
(4rc)2’ 


Ne = t, z(£ = 0 ) = y, 


z(£ = t) = x . 


Here we used the path integral Feynman-Schwinger representation. Thus, the solution of 
(1) has the form 


L 

<p(x,t) = j dt' J dy(Dz) xy K z (t,t')F z (t,t'}g(y,t') , 


where 


K z (t, t') = 6(t — t')exp | —m 2 {t — t') — J ——— dt; J , 

V t’ 

F z {t,t') = exp (g j <p{z(£),£)d{ > 

V v > 

To obtain the system of eqs. for < cp >, < iprj >, < pp >, < ppg >,... we use the 
generating functional 

Z[J] =< F z (t, t')F z (t, t')exp (^J dudtJ(u,t)rj(u,t> 

and apply to it cumulant expansion [6]: 


oo 2 

Z[J] = exp^2 — < 


/ * 


n =i m v . 


9 / <p(z(€),Qd£ + g / + / dudtJ (u, t)rj(u, t) 


> . 


Then 


< F z (t, t')r)(y, t') >= 


5Z[J] 

8 J(y,v) 

n —1 


J =o — 

Fz(t,t ')=1 
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< (gj ¥>(*($)> v(y,t') > J < ^(M') > , 
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< F z (t,t')r](y,t')ri(u,t) >= 


S 2 Z\J\ 


6J(y, t')5J(u, t ) 


n —2 


j=0 

Fz(t,t')=1 


oo ^ 
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< F z (t,t') > 


and so on. In the bilocal approximation (i.e. omitting all the cumulants higher than 
quadratic), one gets: 

< (p(x, t) > = 

t t 

= 9 1 dt> I d y( Dz )^v K z(t,t') f d £ < >< ^(M') > , ( 2 ) 

0 t' 

< ip{x, t)rj(x, t) > = 


= 2 


L 

J ( Dz) xx K z (t,t) < F z (t, t) > +g 2 J dt' J dy(Dz) xy K z (t,t')- 


t t. 

j d £< <p(z(Q,Qv(x,f) > f d £ < ¥>(*(£')> >< F z (t, t') > , (3) 

t' t' 

< (f(x, t)ip(x , t) > = 

t t 

= g 2 J dd J dt 1 J dydy(Dz) X y(Dz) X yK z (t,t , )K z (t,i')- 


0 0 


i L 

J d t< <p(z(€),€)v(yj) >+ j d£' < <p(z(g),g)ri(y,t') 


> 


t t 

f d£ < <p(z(£),£)rj(y,F) > + J d£' < <p(z(£'),£')ri(y,V) > 

1 1> t 1 

t 

■ < F z (t,t')F z (t,t') > +2 f dt' f dy(Dz) xy (Dz) xy K z (t,t')K z (t,t') < F z (t,t')F z (t,t') >, 


o 


( 4 ) 


where 


<F z (t,t')> = =exp [g d£ <tp(z(£),£) > + 


3 



2 t t 

+ !/>/ d('(< AAAAM4('),C) >-< vWM) >< A40.0 >)) , 

t' t' 

< F z (t,lf)Fg(i,V) > = exp | g J cl£ < (p(z(£),£) > + 
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( 6 ) 


Eqs.(2)-(6) is the minimal closed set of eqs. for < ip >, < ipr\ >, < tpp >. 

Note, that < ip(x)ip{x') > mc = hm t-» < p(x, t)ip{x\ t') > v . It means, that 

t / —► oo 

(t—t') fixed 

without only suppositions about the structure of the real vacuum, we obtained eqs. for 
physical correlators. 


3 Perturbative expansion of eqs. (2)-(6) 

Expanding the right hand size of eqs.(2)-(6) in powers of g, we get in the lowest order 

< (p(x, t ) >( 0) = 0 , 


< c p(x,t)rj(x,t) > (0) = 2 f(Dz) xS K z (t,i) = 


%^)+ m2 d-t) 


8ll 2 (t — t) 2 

The last term on the right hand size of (4) can be written as 


J dy < x\KF\y >< y\KF\x > = < x\KFKF\x > , 


and for t — t one obtains for this term 

2 1 


dti{Dz)^K z (ti,Q) < F z (t i,0) > 


o 


where 


F z (t i,0) = exp 


2 ^ 

J <p(z(£),£ + t-d£ + j <^(<z(£),£ + t — tfjdf 


( 7 ) 
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Note, that in the asymptotical regime, when one drops out the dependence on £ in 
^Cz(£)>£)> we liave F z F z • Using (7), 

< lf(x,t)(f(x,t) > ( '° ) = 

dp e ip{x-x)(\ _ g-2(p 2 +m 2 )t\ 


2 1 


J dti J (Dz) xl Kz(t u 0) = j 


(2tt) 


p2 _|_ m 2 


This expression tends to the propagator of a free boson in the limit t —> oo. 
In the next order we have 


= 9 



that in the limit t —» oo corresponds to the lowest order tadpole diagram. 
Keeping < p in the last term on the right hand size of (4), we obtain 

, / /- (1) f 1 , f dp ip{x-x)-(p 2 +m 2 -g<ip> ( - 1 '>')t 1 

< <p(x,t)<p(x,t) > ( >= J dtij J^y e V J = 


[ d P eW{X , 0fe 2 to< V >W-m*)U 

J (27r) 4 p 2 + m 2 — g < p >d) - 


( 8 ) 


Expanding (8) in powers of g, one obtains an infinite set of diagrams. 

Note, that the unstability of vacuum in the p 3 theory leads to divergence of (8) if 
g < p >( 1 i> m 2 , so that the stochastic process has not limiting equilibrium. 


4 Gluodynamics 

We shall start from the Langevin equation 

K = (.D A F,„)“ - (9) 

where F ^ = d\A“ - d r A\ + gf^A^A^, = 9aF jJ, + gf" J " : A l [F^ / , and the sign of 

? 7 “ is changed. Let’s use Schwinger gauge [7] A ll {x,t){x — x 0 )/t = 0, in which A^(x,t) = 
f* 0 dz„a(z,x)F u/i (z,t), where a(z,x) = , %o is an arbitrary point (here and 

later in all the integrals of the type dz v the path of integration is a straight line), and 
introduce the generating functional 

$p{t) = P exp ig cte M ( f dz u a(z, x)F Uft (z, t ) + t)) 

J c J x o 
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, where n p (x,t) = Jo Vn(x, t')dt', C is some fixed closed contour. 
According to (9), 



< $p(t) >= igtr 


du p < $p(t)V p ((3,u,x 0 ,t) >, 


( 10 ) 


where 

V p ((3,u,x 0 ,t) = <h(a;o ,u,t)(D x F Xfl (u,t) + (/3 - l)rj p (u, t))${u, x 0 , t), 

ru 

<h(u, x 0 , t) = P exp ig / A p (z, fjdz^. 

Jxo 

Applying to both sides of (10) cumulant expansion, using the formula [6] 

OO 1 

< e A B >=< e A > (< B > + Y — < A n B »), 

yi nl 

where A and B are two statistically dependent matrixes, and putting (3 = 1, one obtains 
in the bilocal approximation 

rv r u d fy r t d 

try / dz x a(z,y) dx p a(x,u)—< F Xu {z,Xo,t)F pp (x,x 0 ,t) > + dz x a(z,y )/ dt —■ 

Jxo Jxo Ot Jxq Jo Ot 


ru rt Q 

< F Xu (z , x 0 , t)rj p (u, x 0 , t, t')>+ dx p a(x , u) / dt'— < F pp (x, x 0 , t)rj„(y , x 0 , t, t') > + 

Jx o Jo Ot 

rv 

+ / dz x a(z, y) < F Xu (z, x 0 , t)g p {u, x 0 , t,t) > + 

Jxo 

ru 

+ / dx p a(x, u) < F Pfl (x, xq, t)y v (y, x Q , t,t) > + 

Jxo 

+ dt\< r] u {y, x 0 , t, t)r) p (u, x 0 , t, t') > + < r) v {y, x 0 , t , t')g p (u, x 0 , t, t) >) + J dt' dt"■ 

// py c) 

< r)v(y,x 0 ,t,t')r) p (u,x 0 ,t,t") >) = 2 tr( / dz x a(z,y)— < F Xu (z,x 0 ,t)F pp (u,x 0 ,t) > + 

JXQ Ulin 


d_ 

"dt 


rt , Q 

+ / dt — < F pp (u,x 0 ,t)r] w (y,x 0 ,t,t') >), 
Jo ou 0 


( 11 ) 


where 


F Xv (z, x 0 , t) = §(x 0 ,z,t)F Xu (z,t)$(z,x 0l t),ri P (u,xo,t,t') = ®(x 0 ,u,t)- 


■rj tl (u,t')$(u,x 0 ,t). 

Differentiating (10) by (3 , we get 

d / 

tr— < <& 0 (t)n p (u, x 0 , t.) >= tryig j) dz u < ®p(t)n u (z, x 0 , t)V p ((3 , u, x 0 , t) > + (12) 

+ < ^y(t)Vu,(u,x 0 ,t) >), 
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that in bilocal approximation yields 



dz x a(z,y ) < F Xu (z, x 0 , t)n p (u, x 0 , t) > + < n u (y,x 0 ,t)n p (u,x 0 ,t) >) = 


= tr(< $(x 0 ,u,t)(D x F XlJ/ (u,t))$(u,Xo,t)n u (y,Xo,t) > + 



dz x a{z,y ) < F Xu (z,x 0 ,t)- 


■y n(u, x 0 , t) > + < n v {y , x 0 , t)y p (u, x 0 , t ) >). 


Noticing that 


tr(< <3>(. x 0 ,u,t)(D x F Xp ,(u,t))^(u,x 0 ,t)n I/ (y,x 0 ,t ) >= 

d f u 

= o —tr < F XfX (u , x 0 ,t)^?/, x 0 , t) > +ig / dx p - 

OU\ J xo 

•ot(x, u)(tr < F Xfi (u,x 0 ,t)n v (y,x 0 ,t)F X p(x,x 0 ,t ) > -hr < F A/i (w, x 0 ,i)F Ap (x, x 0 , t)- 

• n u (y,x 0 ,t ) >) and using the expression for threelocal path-ordered cumulant -C 
123 >=< 123 > - < 1 >< 23 > - < 12 >< 3 > + < 1 >< 2 >< 3 > [7], one 
obtains 


tr( J dz x a(z,y )J^ < F Xv (z, x 0 , t^K x 0 , t, t') >+^ dt"j^ < y u (y, x 0 ,t')r] p (u, x 0l t,t") > 


d 


- w— < £ 0 , h)^(y, £ 0 , t, t) >) = -tr < y u (y, x 0 , t, t)y p (u, x 0 , t, t') > . (13) 

ou x 

Differentiating (12) by /?, putting (3 = 1 and using the expression for fourlocal path- 
ordered cumulant -C 1234 S>=< 1234 > — -C 123 S> < 4 > — < 1 ><C 234 S> — -C 
12 X 34 > - < 12 >< 3 >< 4 > - < 1 >< 23 X 4 > - < 1 >< 2 >< 34 > - < 
1><2><3><4>, we have in bilocal approximation: 


tr 


0 


t ,, d 

dt — < rj v (y, x 0 , t, t')y p (u , x 0 , t, t) > + < y u (y, x 0 , t, t)y p (u, x 0 , t, t') > - 




ru 

- < 7}„(y, x 0 , t, t')y p (u, x 0 , t, t) >) = g 2 tr( / dx a a(x, u ) < F w (u, x 0 , t)F pa (x , x 0 , t) > • 
j>dz x dt" < i q x {z,XQ,t,l!)'n v {y,x Qi tJ l ) > - £dz x j dt' < F plt (u,x 0 ,t)rix(z,xo,t,lf) > 

f u n 

■ / dx a a(x,u ) <r) v (y,x 0 ,t,t )F pa (x,x 0 ,t) >) (14) 

Eqs. (11),(13) and (14) is the minimal closed set of eqs. for correlators 


< F Xu (z, x 0 , t)F pp (x, x 0 , t) >, < F Xu (z, x 0 , t)y p (u, x 0 , t, t') >, 

< r] u (y,x 0 ,t,t')r] p (u,x 0 ,t,t') > . 

Note, that beyond obtained eqs. (and all other eqs., one is able to get, using 
<&p) there are some additional relations, connecting gauge-invariant correlators. Let 
G w ... Mn (xi,ti, ...,x n ,t„,x 0 ,t) be a product of some number of F pu and (or) rj M , which 
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are given in the points x\,.,, r x n at the moments t±, ...,t n of fictitious time respectively, 
and all the parallel transporters between x$ and each of these points are given at the same 
moment t. Then, using nonabelian Bianchi identities D X F Xp = 0, we have [7]: 

d r x 

tr -— < >= Wtr / dz p a(z,x)(< F Xp (x, x 0 , t')F Xp (z, x 0 , t)G Pl ... Pn > - 

OX\ Jx o 


- < F Xp (x, x 0 , t')G Pl ... Pn F Xp (z, x 0 , t ) >), where x A x 1: ..., x ± x n . 

To check ourselfs, let’s consider eqs. (11), (13) and (14) in the order g°. Looking for 
< A“(t/, t)r] b ^(u, t') > in the form S ab d pu (z, r), where 

z = u - y,r = \t - t'\, d up (z, t) = d flu (z,r),d llu (-z,T) = d pi/ (z,r), one obtains from (13): 
(d pX -§^ + dz ^ dzx - 5 pX d 2 )d Xu {z,r) = -2 5 pu 5(z)5{r) and, thus, 

dpv(k, r ) = —2 d{r)(Tp U e~ k2T + L pu ), where T pu = S pu - L pu = ^F. 

Looking for < A^(y,t)A b fl (u,t) > hr the form 5 ab h pv {zA)^ where 
K p (z,t) = h lxu (z,t),h l _ tu (-z,t) = h pu (z,t), we have from (11): 

(|<WS + - S pX d 2 )h x „{z,t) = -d pu (z,0). 

So, h p „(k,t) = 1 — e ~ 2k *) + 2 tL pu , that is the ordinary photon propagator in the 

stochastic quantization method [1]. 


5 Conclusion 

The main result of the present letter is eqs. (11), (13) and (14). They are explicitly gauge 
invariant and produce correct perturbative results in the lowest order. 

The principle of separation of perturbative and nonpertnrbative contributions in the 
obtained eqs., introduction of quarks and the problem of regularization will be treated in 
the next publications. 

One of us (Yu.S.) would like to thank Yn.M.Makeenko for useful discussions. 
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